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Chimeras in networks with purely local coupling 

Carlo R. LaingQ 

Institute of Natural and Mathematical Sciences, Massey University, 

Private Bag 102-904 NSMC, Auckland, New Zealand, 
phone: +64-9-414 0800 extn. 43512 fax: +64-9-44$ 9790 
(Dated: June 22, 2015) 

Chimera states in spatially extended networks of oscillators have some oscillators synchronised 
while the remainder are asynchronous. These states have primarily been studied in networks with 
nonlocal coupling, and more recently in networks with global coupling. Here we present three 
networks with only local coupling (diffusive, to nearest neighbours) which are numerically found to 
support chimera states. One of the networks is analysed using a self-consistency argument in the 
continuum limit, and this is used to find the boundaries of existence of a chimera state in parameter 
space. 
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Chimera states, in which a symmetric network of iden¬ 
tical oscillators splits into two regions, one of coherent 
oscillators and one of incoherent, have been studied in¬ 
tensively over the past decade |1 3j. Spatial networks on 
which they have been studied include a one-dimensional 
ring 0,0i3La square domain without periodic boundary 
conditions [til Hlj , a torus 1313 and a sphere 0. They 
have also been observed recently in a number of experi¬ 
mental settings 00. Early investigations considered 
networks with nonlocal coupling (i.e. neither all-to-all 
with uniform strength, nor local coupling, via diffusion, 
for example) as chimeras were first reported in nonlocally 
coupled systems 03. Nonlocal coupling was at first 
thought to be essential for the existence of chimeras, how¬ 
ever, more recent results show that chimeras can occur 
in systems with purely global coupling 30. 

Here we consider the opposite limit and address the 
question as to whether chimera states can exist in spatial 
networks with purely local coupling. We present three 
such networks in which this does occur. The idea behind 
the creation of the networks is straightforward and can 
be found in the early papers [ll|, [2(J. Consider a gen¬ 
eral reaction-diffusion equation on a one-dimension spa¬ 
tial domain Ul with only local interactions via diffusion 
in one variable: 

=/(«)+« (1) 

dv d 2 v 

e m= 9(u) - v + d. ? <2> 

Setting e = 0 in (|2]) one has 

0 -1?)"= 9(u) (3 > 

and if h(x ) is the Green’s function associated with 



on Q then we can write (j3j) as 


v(x) = h(x- y)g(u(y)) dy (4) 

Jn 

and substituting this into 0 we obtain a closed nonlocal 
equation for u. We will implement the network analogue 
of (P)-® but with e small and nonzero in the expectation 
that the behaviour of interest when e = 0 persists for 
e > 0. 

The first model we consider consists of N oscillators, 
equally-spaced on a domain of length L, with periodic 
boundary conditions. The state of oscillator j is de¬ 
scribed by two variables: 0j G [0, 27r) and Zj G C. (A 
complex variable is used to simplify presentation; we 
could equally well use two real variables.) The governing 
equations are 

1th =Uj- Re( Zj e- i09 ) (5) 

c dz i z- I *J+i ~ 2 *i + *i-i (6) 

dt~ 3+ (Ax) 2 W 

for j = 1,2... TV, where Ax = L/N and A, /3 and 
e are all constants. The ujj are randomly chosen 
from a Lorentzian distribution with half-width-at-half- 
maximum a centred at ojq, namely 


sM = 


CF / 7T 

(cj - uj 0 ) 2 + cr 2 


( 7 ) 


An example of the system’s dynamics are shown in Fig.P 
The domain clearly splits into two regions, one showing 
coherent behaviour of the phases and the other, incoher¬ 
ent. This behaviour has been replicated in networks of 
up to N = 1000, so is not a small-AT effect. 

To understand the relationship between ©-© and 
previously studied models we set e = 0 in (j6]). If Zj 
is the jdh entry of the vector z G C N and similarly for 0j 
we can write © as 


(/ - D)z = Ae i(e+0) 
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To analyse the chimera seen in (j5j)-(j6l) we use a self- 
consistency argument similar to that in [1,(2,1 5| . We first 
move to a rotating coordinate frame, letting 0 j = Oj — fit 
and y = Zje~ lQt , where fl is to be determined, and then 
take the limit N —>■ oo, to obtain 


— = w - Re [ye l<t> ) - 9 ( 11 ) 

-y+jJL- ie9y (12) 

We now search for solutions of CCED-CED for which y is 
stationary, i.e. just a function of space. We let such a 
solution be y = R(x)e l@ ( x \ Since y is constant we can 
use m to determine the dynamics of 0 for any y and 
uj\ if \R\ > |cj — fl\ then 0 will tend to a stable fixed 
point of (HU whereas if \R\ < \oj — f2| then 0 will drift 
monotonically. To obtain a stationary solution of m 
we replace e^ by its expected value, calculated using 
the density of 0 , which is inversely proportional to its 
velocity (given by <m»- So (keeping mind that u) is 
random variable) we need to solve 


0 = Ae ip 


/ oo p 2 ir 

/ ' 

-oo J 0 


uj)g{uj)d(j) duj — y - 


&V 

dx 2 


where the density of 0 given uj is 


iefly 

(13) 


Figure 1: (Color online) Chimera solution of the system ©- 
©. (a): sin0 7 : (b): sin (arg (zh)): (c): \zh\. Parameters: 

LJo = l,a = 0.01, e = 0.2, A = 1.5, L = 2tv, N = 100, p = 0.1. 


where I is the N x N identity matrix and D is the matrix 
representation of the classical second difference operator 
on N points with periodic boundary conditions. Defining 
G = (I — D)~ 1 we have 


N 

Zj m A'£ l G jk e i V* + n (9) 

k=l 

where Gjk is the jkth element of G, and substituting © 
into © we obtain 



N 

A ^2 Gjk cos (9j -0 k - (3) 

k=l 


( 10 ) 


Since I — D is circulant, so is G, i.e. G^ is a function of 
only \i — j\, and (flOl) is thus of the same form as studied 
by a number of others HiiB- An important property 
of ([TO]) is that it is invariant with respect to a uniform 
shift of all phases: 0j Oj + 7 for all j where 7 is some 
constant. This implies that (flOl) can be studied, without 
loss of generality, in a rotating coordinate frame where 
cjo = 0 ? he. the actual value of ujq in m is irrelevant. 
This is not the case for (|5])-(|6]) when e 7 ^ 0 (although 
© is invariant under the simultaneous shift: 0j 0j -\- 7 
and Zj 1 y Zje 11 for all j ). 


p{<t>\u) 


vV - o ) 2 - W 

2n\uj — fl — R cos (@ — 0)| 


(14) 


and g(uj) is given by ©. Evaluating the integrals in (fl3l) 
we obtain 


^ e i(©+/3) - 

T 


cjq + icr — D — \/ (w 0 + i(J — fl ) 2 — R? 


1 + iefl 


d 2 \ 
dx 2 ) 


Re ^ = 0 (15) 


We determine 0 and fl by simultaneously solving (fl5l) 
and the scalar equation 0 ( 0 ) = 0 ; this equation amounts 
to choosing the origin of the rotating coordinate frame. 
Following solutions of m as co>o and e are varied we 
find that two solutions are destroyed in a saddle-node 
bifurcation on the curve shown in Fig. [2l Although 
our self-consistency argument gives no information about 
the stability of solutions (unlike the continuum theory 
in [4, L 7 , 3 ) quasistatic sweeps through parameter space 
indicate that the curve in Fig. [2] does indeed mark the 
boundary of stable chimeras in the system ©-©. if e is 
increased past the boundary in Fig. [2] when ujq is to the 
left of the cusp (at ujq « 1.04) the system ([S])-© moves 
to the almost synchronous state, whereas to the right of 
the cusp the system moves to a spatially-disordered state, 
and the almost synchronous state seems unstable here. 

As a second example we consider a network of Stuart- 
Landau oscillators, each of which can be thought of as 
the normal form of a supercritical Hopf bifurcation, with 
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Figure 2: Saddle-node bifurcation of solutions of Hu. 
Chimera solutions of ©-© are stable below the curve. Pa¬ 
rameters: A — 1.5, L — 2i r, P = 0.1, a — 0.01. 


purely local coupling through a second complex vari¬ 
able. Using Stuart-Landau oscillators as opposed to the 
phase oscillators above introduces an amplitude variable 
to the oscillator dynamics. As above we have N oscilla¬ 
tors equally-spaced on a domain of length 1 with periodic 
boundary conditions. The equations are 


— (1 + ivo)Aj - (1 + ib)\Aj\ 2 Aj 

+ ^0- + ~ Aj) 


dZ j 
dt 


— Aj Zj 


Zj+i - ‘ZZj + Zj -1 

16(Ax) 2 


(16) 

(17) 


for j = 1,2 ...N where Aj,Zj G C and uJo,a,b,K and 
e are real parameters and Ax = 1/N. Note that the 
oscillators are identical. A chimera state for this system 
is shown in Fig. [3] Note that setting e = 0 in (H3 and 
then taking the limit N oo we obtain the nonlocally 
coupled complex Ginzburg-Landau equation for just the 
variable A , as studied by [lj]. Then assuming that K is 
small one finds a scale separation between the amplitude 
and phase dynamics of A and upon setting \A\ = 1 the 
phase dynamics can be written in a nonlocally coupled 
form p], |s]. 

As a third model we consider a heterogeneous network 
of oscillators, each described by an angular variable and 
a real variable. The angular variables have the form of 
Winfree oscillators [24l-[26j. The model is 


m j 

dt 


= ujj + K J Q(0j)uj 


duj _ N 
e ~dt ~ _ U J 


U j _|_ \ 2 Uj | U j — 2 . 

(. Ax ) 2 


(18) 

(19) 


for j = 1, 2 ... N where Q(6) = sin ft — sin (0j + /3), ft, f3 
and e are parameters, P n {6) = a n ( 1 + cos^) n where 
n > 1 is an integer and a n = 2 n (n!) 2 / (2n)\ (so that 



Time 



Time 


- 0.14 


0 - 0.145 

0 

=3 

cr 
o 


- 0.15 


- 0 . 155 . 


(C) 

. "X 

V 

% 

. m* . 






0 


50 


100 

Index 


150 


200 


Figure 3: (Color online) Chimera solution of the system (1161) - 
(H7I) . (a): \Aj\; (b): | Zj\; (c): average rotation frequency of 
the Aj over a simulation of duration 2000 time units. Param¬ 
eters: cjq = 0, e = 0.01, a — —1, b = 1, K = 0.1, N m 200. 


/ 0 27r P n (0)d0 = 27r) and Ax = L/N. The ujj are ran¬ 
domly chosen from a normal distribution with mean co>o 
and standard deviation a. Q{6) is the phase response 
curve of the oscillator and can be measured experimen¬ 
tally for a neuron, for example [27j . 

A chimera state for m-m i s shown in Fig. 2J Set¬ 
ting e = 0 in m and solving for the Uj one would ob¬ 
tain a nonlocally coupled network of Winfree oscillators. 
Chimeras have been found in a network of two popula¬ 
tions of Winfree oscillators [24|, [28] but a truly nonlocally 
coupled network has apparently not yet been studied. 

We have presented three one-dimensional networks, 
where each node is described by one variable which has 
a phase associated with it and a second variable which 
is coupled in a diffusive fashion to just its two nearest 
neighbours. All networks have the same structure and 
show chimera states over some range of parameters. All 
have a small parameter (e) which controls the time scale 
of the diffusing variable, so can be thought of as slow-fast 
systems [29| . 

We have not given any stability analysis of the models 
presented here, only a self-consistency argument for the 
first model. Chimeras in systems of the form ([TUI) have 
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been studied by passing to the continuum limit (N —>> oo) 
and analysing the resulting continuity equation using the 
Ott/Antonsen ansatz [1, |8j, 30, l3lj | . However, it does 

not seem that such an approach could be used to study 
the models presented here due to the dynamics of the 
extra variables. 
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Figure 4: (Color online) Chimera solution of the system (H8|) - 
G9- (a): sin6b; (b): Uj\ (c): average rotation frequency of 
the 6j over a simulation of duration 2000 time units. Pa¬ 
rameters: ujq — 0.3, cr — 0.001, n s ® 4, L = 4, k «= 0.4, /3 m 
7t/2 — 0.2, e = 0.1, N = 100. 


Regarding experimental implementation, note that the 
nonlocal coupling in the experiments reported in [16 1 2 3 4 5 6 7 8 9 10 11 12 , 18] 
was implemented by computer, i.e. the experiments were 
hybrid physical/computer. The models presented here — 
while being caricatures of physical systems — have only 
local, nearest-neighbour diffusive-like coupling. Since dif¬ 
fusion is ubiquitous in spatially extended systems of re¬ 
acting chemicals, the most natural system in which to 
implement networks of the form discussed here (without 
a computer) may be in arrays of microsopic chemical os¬ 
cillators PQ. 
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